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Abstract — In this paper we classify all extremal and s-extremal 
binary self-dual codes of length 38. There are exactly 2744 
extremal [38, 19, 8] self-dual codes, two s-extremal [38, 19, 6] 
codes, and 1730 s-extremal [38, 19, 8] codes. We obtain our 
results from the use of a recursive algorithm used in the recent 
classification of all extremal self -dual codes of length 36, and from 
a generalization of this recursive algorithm for the shadow. The 
classification of s-extremal [38, 19, 6] codes permits to achieve the 
classification of all s-extremal codes with d = 6. 

Keywords: Classification, recursive construction, ex- 
tremal, self-dual codes, s-extremal, shadow 



I. Introduction 

Self-dual codes are one of the most interesting classes of 
linear codes. They have close connections with group theory, 
lattice theory, design theory, and modular forms. It is well 
known that self-dual codes are asymptotically good 11221 . There 
has been an active research on the classification of self-dual 
codes over finite fields and over rings in general (see [25], [23 1 
for details). In particular, the classification of binary self-dual 
codes was started by Pless ll24l and has been actively studied 
by many authors (see lfl9ll for a survey of optimal self-dual 
codes over small alphabets). 

Recently, using a recursive method, Aguilar and Gaborit 
classified all 41 extremal [36, 18, 8] binary self-dual codes. 
These results were pushed further by Harada and Mune- 
masa 1171 who, besides the 41 extremal codes of Qj, also 
give a complete classification of all self-dual codes of length 
36. 

A natural question is hence to consider the case of length 
38. A simple computation on the mass formula shows that 
there are at least 13,644,433 inequivalent binary self-dual 
[38, 19] codes fl7l . It is hence natural to consider the case of 
special subclasses of self-dual codes. The most interesting such 
subclass is the class of exttemal codes. Given the classification 
of all [36, 18, 6] self-dual codes of [ 17 1, we apply an optimized 
recursive algorithm as in JT| to derive the classification of all 
2744 extremal self-dual [38, 19, 8] codes. 
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Another subclass of interesting self-dual codes with com- 
binatorial properties is the class of s-extremal codes: these 
codes are self-dual codes whose weight enumerator is uniquely 
determined, depending on the condition on a high weight 
of the shadow. The notion of codes (and lattices) with long 
shadows was first developed by Elkies ifTTI . This notion was 
generalized by Bachoc and Gaborit in [2J who introduced the 
notion of s-extremal codes. These codes exist depending on 
conditions on their length and their minimum distance. The 
classification of s-exttemal codes with d = 4 was done by 
Elkies. The case of d = 6 was mainly considered in [2], but 
two lengths remained to be classified. One is length 36, which 
was classified in (TJ, and the other is length 38, which is 
what we classify in this paper. Our classification is based on 
a generalization of the subtraction algorithm in the case of 
the shadow. It permits us to use the recursive algorithm by 
showing that in certain cases for n even, the subtraction of 
(11) from a [2n + 2, n + 1, d + 2] self-dual code with shadow 
weight s + 1 leads to a [2n, n, d] self-dual code with shadow 
weight s. This result is interesting in itself. 

The paper is organized as follows: Section II gives pre- 
liminaries and background for self-dual codes, Section III 
compares the different method to extend a self-dual code in 
a purpose of classification. In Section IV we show that there 
are exactly 2744 extremal [38, 19, 8] binary self-dual codes. 
In Section V we prove that there are only two s-extremal 
[38,19,6] codes and 1730 s-extremal [38,19,8] codes. The 
last section describes the covering radii of self-dual codes of 
length 38. 

II. Preliminaries 

We refer to ll20l for basic definitions and results related 
to self-dual codes. All codes in this paper are binary. A 
linear [n, k] code C of length n is a fc-dimensional subspace 
of GF(2) n . An element of C is called a codeword. The 
(Hamming) weight wt(x) of a vector x = (sci, . . . , x n ) is the 
number of non-zero coordinates in it. The minimum distance 
(or minimum weight) d(C) of C is d(C) := min{wt(x) | x G 
C,x ^ 0}. The Euclidean inner product of x — (x\ ,...,!„) 
and y = (y 1 , . . . , y n ) in GF(2) n is x • y = £™ =1 x^. The 
dual of C, denoted by C 1 - is the set of vectors orthogonal 
to every codeword of C under the Euclidean inner product. If 
C = C , C is called self-dual. A self-dual code is called Type 
II (or doubly-even) if every codeword has weight divisible by 
4, and Type I (or singly-even) if there exists a codeword whose 
weight is congruent to 2 (mod 4). 

Two codes over GF(2) are said to be equivalent if they 
differ only by a permutation of the coordinates. Let C be a 
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binary self-dual code of length n and minimum distance d(C). 
Then d(C) satisfies the following (see |25l ). 



d(C) < 



f 4 


n 




.24. 




n 




.24. 



if n ^ 22 
if n = 22 



(mod 24), 
(mod 24). 



A self-dual code meeting one of the above bounds is called 
extremal. A code is called optimal if it has the highest possible 
minimum distance for its length and dimension. 

By the well known Gleason's theorem, the weight enumer- 
ator Wc{x, y) of a Type I code can be written as follows (for 
rational coefficients Cj): 



W c (x,y) 



[n/8] 

5>-' 

i=0 



An important notion associated to a Type I code is the 
shadow S of a code C, defined by S = Cq\C, where C 
is the doubly-even subcode of C. In J9), Conway and Sloane 
show that for a weight enumerator Wc(x, y) given above, the 
weight enumerator Ws of S satisfies 



[n/8] 



4\2i 



i=0 

This notion of shadow permits to give more information on 
potential weight enumerators of self-dual codes, and is also 
used to define s-extremal codes (see j2| or Sec. IVT ). 

The main tool to classify self-dual codes is based on the 
so-called mass formula. It is known from |24| that self-dual 
binary codes (Type I or Type II) of length n satisfy a formula 
(a mass formula): 



N{n) 



x - n! 

*f \Aut{C 3 )\-- 



where the sum is made over all inequivalent self-dual codes 
(Type I or Type II) of length n, \Aut(C)\ denotes the order 
of the automorphism group of a code C, and N(n) is the 
number of Type I or Type II codes. In particular, for Type I 
codes, N(n) = UfJi^ 1 + l ) an d f »r Type II codes N(n) = 

nto 2 (2 i + i). 

Therefore, for n = 38, 
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^(38) = Y[(2 l 



D = 



E 



38! 



\Aut(C 3 )\- 



Hence, 

13644432.20346 



< 



< 



38! , 



\Aut{Cj)\ 

#(all inequivalent self-dual codes) 



Moreover, as there is no mass formula for extremal self-dual 
codes, it might be also difficult to classify all extremal binary 
[38, 19, 8] codes. However, using the recursive construction [l] 
which was used in classifying all extremal binary [36, 18,8] 
codes, we are successful in classifying all extremal binary 
[38, 19, 8] codes. 



A very interesting tool for self-dual codes is the subtraction 
procedure of (1 1) on two coordinates of a code. This procedure 
permits to construct a [2n, n, d! > d] self-dual code from a 
[2n+2, n+1, d+2] self-dual code. It works as follows: suppose 
one starts from a [2n + 2, d + 2] self-dual code C for 
d > 2. Let i and j be two different coordinates of the columns 
of C. Since d+2 > 3 and C is self-dual, any two columns of C 
are independent (if not, there should be a codeword of weight 2 
in C, a contradiction). This implies that the coordinates of the 
two columns of the codewords of C contain (00), (10), (01) 
and (00). For the subtraction procedure of (11) on columns i 
and j, one first keeps all codewords which are either (00) or 
(11) on columns i and j, and then deletes columns i and j for 
these codewords. Let C' be the obtained code. Since d+2 > 2 
and by an argument similar to the shortening of a code, the 
dimension of C is n. Moreover since the scalar product of 
any two codewords of C is 0, the scalar product of any two 
codewords of C is also 0. Now as the minimum distance of 
C is d + 2, the minimum distance d! of C is either d either 
d + 2 (depending on the fact that columns i and j intersect 
or not with codewords of C of weight d + 2). Overall C is a 
[2n, n, d 1 > d] self-dual code. 



III. Construction methods 

There exist several methods to construct self-dual codes of 
length n+2 from self-dual codes of length n. In this section we 
recall these methods; the recursive construction, the building- 
up construction and the Harada-Munemasa construction. We 
eventually compare them. 

A. The recursive construction 

In ITJ, Aguilar and Gaborit give a recursive construction 
of binary self-dual codes. This algorithm can be seen as the 
reverse operation of the subtraction procedure of (11) given 
above. We recall that a subtraction procedure produces a self- 
dual [2n,n,d' > d] code C from a self-dual [2n + 2,ri + 

1, d+2] code C. The recursive algorithm starts from a self-dual 
[2n, n, d] code C and constructs (up to permutation) all self- 
dual [2n+2, n+1, d+2] codes which by subtraction of (11) on 
certain two columns give the code C . The idea of the recursive 
algorithm is very simple and consists of extending the code 
C with 11 for all codewords of weight d, then constructing 
all possibilities with (00) or (11) for a basis of remaining 
codewords, and eventually checking for addition of a vector 
strictly contained in the shadow of the extended code. This 
approach is very useful in classifying extremal self-dual [2n + 

2, n + 1, d + 2] codes because it is sufficient to know (up to 
permutation) a classification of [2n, n, > d] self-dual codes. 
Indeed, any [2n+2, n+1, d+2] code gives a [2n, n, d] code by 
subtraction of (11) on adequate columns, conversely applying 
the 'reverse subtraction' procedure to the set of all [2n, n, d] 
codes (up to permutation) permits to construct a set of codes 
which contains (up to permutation) all [2n + 2, n+1, d + 2] 
codes. 

We now recall the recursive algorithm (with a correction of 
n — k in Step 2) from [1| into f — fc: 



3 



Recursive algorithm 



Input: S n , the set of [n, %,d\ self-dual codes up to permuta- 
tion 

Output: The set of [n + 2, § + 1, d + 2] self-dual codes 
For each code C n of S n do: 

1) List all the words of weight d and construct the subcode Cd of 
dimension k generated by these words. Construct a generator 
matrix Gd of Cd composed only with words of weight d. 

2) Let E be a code of dimension — — k with generator matrix 
Ge such that C n — Cd + E, constructs the extended codes C 
with generator matrices: 



1 

Ol 



1 

Ol 



Gd 



G E 



(1) 



3) 



such that cij G {0, 1}, (1 < i < ~ - k). 
Complete all the previous codes C by nonzero elements of 
C ± /C in order to obtain a self-dual code D and check for 
codes with minimum distance d + 2. For codes with weight 
d + 2 check for the equivalence with already obtained self-dual 



[n + 2, | + 1, d + 2] codes. 



The main result of H) is the following: 

Theorem 1: Applying the previous recursive algorithm to 
the set of all inequivalent (up to permutation) binary self- 
dual [n,n/2,d] codes permits to find all inequivalent self-dual 
binary [n + 2, n/2 + 1, d+ 2] codes. 



B. The building-up construction 

There are other constructions generating self-dual codes of 
length n + 2 from self-dual codes of length n. In particu- 
lar, we compare the above construction with two construc- 
tions; the building-up construction |21| by Kim, and Harada- 
Munemasa's construction ifTTl since both constructions gen- 
erate all self-dual codes of length n + 2 from the set of all 
self-dual codes of length n. 

Theorem 2: ([21, building-up]) Let Go = (rj) be a 
generator matrix (may not be in standard form) of a self-dual 
code Co over GF{2) of length n, where is a row of Go 
for 1 < i < n/2. Let x be a vector in GF{2) n with an odd 
weight. Define yi :— x • Yi for 1 < i < n/2, where • denotes 
the usual inner product. Then the following matrix 





1 





X 










G = 














Go 




Vn/2 


Vn/2 





(2) 



generates a self-dual code G over GF(2) of length n + 2. 
The converse of the building-up construction holds as fol- 
lows. 



Theorem 3: (||2TJ) Any self-dual code G over GF{2) of 
length n with minimum weight d > 2 is obtained from some 
self-dual code Go of length n — 2 (up to equivalence) by the 
construction in Theorem [2] 

The recursive construction is a special case of the building- 
up construction. The reason is as follows. 

We show that the matrix in the form (TJ together with a 
representative in C /G whose weight is > 2 can be written 
in the form (01 up to permutation equivalence. Suppose we 
are given the matrix in the form ([TJ above and let G be the 
code generated by this matrix. Then there are four cosets of 
G in G x ; that is, G, zi + G, z 2 + G, and z x + z 2 + G for 
some nonzeroes zi,z 2 G GF(2) ,l+2 . We may assume that 
zi = (1, 1, 0, 0, • • ■ ,0) since zi is nonzero and orthogonal to 
C. Then the minimum weight of G U (zi + G) is 2, which is 
excluded. Hence by permuting the first two columns of z 2 if 
needed, we may put z 2 = (1,0 | x) where x G GF(2) n . As 
G U (z 2 + G) is designed to be self-dual, z 2 is orthogonal to 
itself; hence x is odd. Then as z 2 • (1, 1 | r^) = 0, where is 
a row of Gd in the form (fl]i for 1 < i < k, we have x • = 1. 
Thus by letting := x • = 1 for 1 < i < k, we obtain 
the matrix of the form (|2). This implies that the recursive 
construction is a special case of the building-up construction. 



C. The Harada-Munemasa construction 

In what follows, we recall Harada-Munemasa's construc- 
tion [17|. We note that this is a binary version of Huffman's 
construction ifTSIl for Hermitian self-dual codes over GF(4). 

Let Gi be a generator matrix of a self-dual [n, n/2, d] code 
C\. Then the matrix 



G 2 := 



a i 



n/2-l 



a i 



i/2-l 



(3) 



where a% G GF(2) for (1 < i < n/2 — 1), generates a self- 
orthogonal [n + 2, n/2] code G 2 . The matrix of the form (O is 
a general form of (Q~|i in the recursive construction. In order to 
reduce the possibilities of a,'s, they IfTTl consider the orbits of 
the vector a T :— {a±, ■ ■ ■ , a n /2~i) T under a certain subgroup 
of GL(n/2 — 1,2) to get equivalent self-dual codes of length 
n + 2. After reducing the possibilities, as in the recursive 
construction, add to G 2 a coset z 2 + G 2 from Gj 1 / G 2 whose 
weight is > 2 to get a self-dual [n + 2, n/2 + 1, d' > 2] 
code. Unlike the recursive construction, Harada-Munemasa's 
construction does not necessarily give self-dual [n+2, n/2+1] 
codes with minimum weight d' = d+2. 

D. Comparison of the different methods 

The recursive construction is specially interesting when one 
wants to classify extremal codes since it permits to obtain a 
partial classification for a given minimum distance while other 
constructions do need to start from a whole classification. 

More precisely, the recursive construction is more efficient 
than the building-up construction in generating many self- 
dual codes with higher minimum weight. This is because the 
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recursive construction checks a relatively small number of pos- 
sibilities of a^s in Step 2), whose complexity is 2"/ 2 ~ fe , where 
k > 1 depends on the given code. From our experimental 
results, the dimensions k of subcodes of the 58671 [36, 18, 6] 
codes generated by linearly independent vectors of weight 6 
lie between 2 and 18. We give the possible values of k and 
the number num of their subcodes in Table Q] 

TABLE I 

Number of self-dual [36, 18, 6] codes whose subcode generated 
by codewords of weight 6 has dimension k 



dim k 


num 


dim k 


num 


dim k 


num 


2 


148 


8 


4615 


14 


8170 


3 


5 


9 


911 


15 


5311 


4 


666 


10 


7165 


16 


6290 


5 


45 


11 


2299 


17 


4492 


6 


2165 


12 


8411 


18 


3615 


7 


263 


13 


4100 







We see from our table that there are much more subcodes of 
large dimension than those of small dimension and this clearly 
shows the efficiency of our recursive algorithm. 

On the other hand, the building-up construction [21 1 needs 
2" _1 possibilities for the choice of odd vectors x, generating 
all self-dual codes with various minimum distances. This 
complexity can be reduced to 2"/ 2 as remarked in fl3l . which 
is still higher than that of the recursive construction. 

As described above, Harada-Munemasa's construction is 
effective if the given code has a large automorphism group in 
order to reduce the complexity of checking the equivalence. 
For example, if n = 36, then 41019 (respectively 11242) 
out of the 58671 self-dual [36, 18, 6] codes ___ have the 
automorphism group order 1 (respectively 2). Thus Harada- 
Munemasa's construction usually requires 2 19 or 2 18 possi- 
bilities to generate self-dual codes of length 38 with various 
minimum distances, given a [36, 18, 6] self-dual code. 

Overall, we conclude that when we classify binary self-dual 
[38, 19, 8] codes, the recursive algorithm is much faster than 
the other two constructions. 

IV. Classification of the [38, 19, 8] self-dual codes 

A. Construction of all [38,19,8] self-dual codes 

There are two possible weight enumerators W\ , W 2 and 
shadow weight enumerators Si , 52 for an extremal self-dual 
[38,19,8] code __. 



Wi = 1 + 171/ + 1862y 10 + ■ • • 

51 = lUy 7 + 90447/ 11 + 118446i/ 15 
W 2 = 1 + 203y 8 + 1702y 10 + • • • 

5 2 = y 3 + 106y 7 + 9072?/ 11 + 118390y 
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(4) 
(5) 
(6) 
(7) 



In __ two self-dual [38,19,8] codes with W\, denoted 
by i?3 and D±, were given, where |Aut(i?3)| = 1 and 
|Aut(L> 4 ) = 342. In [16) one self-dual [38,19,8] code C 38 
with W 2 was given with Aut(C3s) = 1. Then Harada [15] 
gave 40 self-dual [38, 19, 8] codes with Wi and W 2 and 
automorphism group orders 1,2,4,8. Later, Kim [21] con- 
structed 325 self-dual [38, 19, 8] codes with Wi and W 2 and 



automorphism group orders 1,2,3. Hence there are at least 368 
inequivalent self-dual [38, 19, 8] codes. We show that there are 
exactly 2744 inequivalent self-dual [38, 19,8] codes. 

Starting from the 58671 [36, 18,6] codes of [17|, we apply 
the recursive algorithm of Section IIII-AI The more expen- 
sive part of the algorithm is the inequivalence testing of 
the differently constructed codes. In order to optimize the 
computation we separated the 58671 [36,18,6] codes into 
sets <S36,i of 1000 codes. To each set, we apply the recursive 
algorithm to obtain a list S^ : i of inequivalent [38, 19, 8] codes 
derived from the set 636,;. Each set 638,1 contains a number 
of inequivalent codes. Then we compared all the Sss,i sets 
to eventually obtained a list of all inequivalent [38, 19, 8] 
self-dual codes. This method permits to avoid many costly 
inequivalence comparisons between codes, since separating the 
whole list of [36, 18,6] codes permits to avoid inequivalence 
testing as the S^s.i list starts from an empty list. 

The whole process took about three weeks on a CPU 
2.53GHz computer. 

Now we obtain our main theorem below. 

Theorem 4: There are exactly 2744 inequivalent extremal 
self-dual [38, 19,8] codes. 

In Table HU we describe all extremal self-dual [38, 19, 8] 
codes with respect to their orders, where |Aut(C)| and num 
stand for the order of automorphism group and the number of 
codes respectively. 

TABLE II 

Number of extremal self-dual [38, 19, 8] codes with respect to 
their orders 



|Aut(C)| 


num 


|Aut(C)| 


num 


[Aut(C)| 


num 


1 


2253 


9 


1 


36 


1 


2 


322 


12 


8 


144 


1 


3 


36 


14 


1 


168 


2 


4 


68 


18 


1 


216 


1 


6 


17 


21 


1 


342 


1 


8 


15 


24 


14 


504 


1 



As mentioned above, the previously known self-dual 
[38, 19, 8] codes have automorphism group orders 1, 2, 3, 4, 8, 
and 342. Hence we list several new self-dual [38, 19,8] codes 
Cgg with different automorphism group orders |Aut(C3 8 )| = 
i = 6,9,12,14,18,21,24,36,144,168,216,504 in 
Appendix. To save space, we only give one code for 
each order. We also list C|| 2 which is equivalent to 
the double-circulant code D4 in @. The list of all 
extremal self-dual [38, 19, 8] codes can be obtained at 



http://www.unilim.fr/pages_perso/philippe.gaborit/SD/GF2/GF2I.htm 



B. An up-to-date table of the number of classified optimal 
self-dual codes 

In the following we give an up-to-date table of the classifi- 
cation of optimal Type I self-dual codes, where being optimal 
means that this is the best possible minimum distance among 
self-dual codes of a given length. These codes may not be 
extremal in the classical sense. For instance, an extremal self- 
dual code of length 34 will have minimum distance 8 if exists, 
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TABLE III 

Number of optimal Type I and Type II codes 



n 


d 


num 


n 


d 


num 


2 


2 




22 


6 


1 


4 


2 


j 


24 






6 


2 




26 


6 


1 


8 


4 




28 


6 


3 


10 


2 




30 


6 


13 


12 


4 




32 


8 


8 


14 


4 




34 


6 


938 


16 


4 


3 


36 


8 


41 


18 


4 


2 


38 


8 


2744 


20 


4 


7 









but it is known that such a code cannot exist and the optimal 
minimum distance is 6. The highest length (up to now) for 
which Type I optimal codes are classified is length 38, which 
is done in this paper for the first time. Notice that it is length 
48 for Type II codes. Complete references for the self-dual 
codes can be found for instance in |fl9l and ll23l . except for 
length 38. 

V. Classification of s-extremal codes 

In this section, we classify s-extremal codes of length 38 
and d = 8 together with s-extremal codes of length 38 and 
d = 6 . 

A. s-extremal codes 

The notion of s-extremal codes was introduced by Bachoc 
and Gaborit in [2|. This type of codes is related to the notion 
of self-dual codes with long shadows introduced by Elkies in 
ifTTl . We recall the definition of s-extremal codes from 13. 

Let C be a Type I self-dual binary code of length n. We 
denote by Co the doubly-even subcode of C. We denote by x 
an element of C\Cq. The shadow S is defined by S = C 1 \C, 
we denote by y an element of S\C. We have Cg 1 = C0UC1 U 
C2UC3 for Ci = y + C ,C 2 = x + C and C 3 = x + y + C . 
Then it is well known that C — C U C 2 and S = C\ U C 3 . 
We have moreover the following three facts [9|: 

1) for any y e S, weight(y) = S (mod 4) 

2) for any y E S and x e C 2 : x ■ y = 1, 

3) for any y E S and z £ Co ■ x ■ z = 0. 

We denote the weight enumerators of C and S by Wc and 
Ws, respectively. From [9|, there exist Co, . . . , C[ n /S\ € R such 
that: 

Wc(x, y) = El^o 81 c % {x 2 + 2/ 2 )i- 4 '{zV(z 2 - y 2 fY 
W s (x,y) = El^f 1 Cl (-l)*2f- 6 *(xy)f- 4l (z 4 - y 4 ) 2 * 

(8) 

Let d be the minimum weight of C and s the minimum 
weight of its shadow. 

Theorem 5: (2) Let C be a Type I self-dual binary code of 
length n with minimum weight d, and let S be its shadow 
with minimum weight s. Then, 2d + s < 4+ ||, unless n = 22 
mod 24 and d = 4[ra/24] + 6, in which case 2d + s = 8 + §. 

A Type I code whose parameters (d, s) satisfy the equality 
in the previous bounds is called s-extremal. In that case, the 
polynomials Wc and Ws are uniquely determined. 



A bound for n when the minimum weight d of an s-extremal 
code is divisible by 4 has been given in lfl2l and in fl4l . and 
a bound has also been given for d — 6 [2, Theorem 4.1] and 
d = 2 (mod 4) with d > 6 Q3- 

Theorem 6: (112, lH4l ") 

Let C be an s-extremal code with parameters (s,d) of length 
n. If d = (mod 4), then ra < 6d - 2. 

Theorem 7: ( lfT4l l 

Let C be an s-extremal code with parameters (s,d) of length 
n. If d > 6 and d = 2 (mod 4), then n< 21d - 82. 

Before proving our classification of s-extremal codes of 
length 38, we prove a result which permits in certain cases to 
relate the weight of the shadow of a code C with the weight 
of the shadow of a subtracted code by (11): 

Theorem 8: If C is a [An + 2, 2n + 1, d + 2] self-dual code 
with d = (mod 4), d ^ and shadow weight s > 3, then 
there exist two coordinates of C on which the subtraction of 
(11) gives a self-dual [An, 2n,d] code C with shadow weight 

8-1. 

Proof. Our proof is based on the existence of the following 
four vectors x, y, z, and t such that: 

1) y = (y'10), y e S of weight s 

2) x = (x'll), x e C 2 of weight d + 2 

3) z = (z'll), z S C* 

4) t = (t'10), t e C 

Let y G S of weight s and x 6 C2 of weight d+2. We have 
x-y = 1, that is, x and y meet in an odd number of positions. 
Then Xi = 1 = y,: for some i. As the weight of x is even, 
there is a j such that Xj = 1 and yj =0. Up to permutation, 
we may assume that x = (x'll) and y = (y'10). Now it 
remains to show that there exist z and t given above. To do 
this, note that Cq = C U S. Hence the minimum distance of 
Cq 1 = min{d + 2, s} > 3. Hence every two columns of a 
generator matrix of Co are linearly independent. (This means 
that C has strength 2. See EU] p. 435] for the term.) Thus 
in each set of two columns of Co each binary 2-tuple occurs 
the same number |Co|/4 of times. Therefore there exist z = 
(z'll) e C and t = (t'10) £ C . 

Since the coordinates of z and t are respectively (11) and 
(10) on the last two positions, there exists a doubly-even code 
Cq of dimension 2n — 2 such that the doubly-even subcode 
Co of C can be written: 

/ z' 1 1 \ 
t' 1 


eg i ! 
V 0/ 

Now if one subtracts (11) on the two last columns of C one 
obtains a code C', such that its doubly-even subcode C has 
dimension 2n — 1, (2n — 2 vectors of Cq and the vector x' - 
which cannot be null since d + 2 =^ 2), the subcode Cq can 
be written as: 
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Overall a generator matrix of C can be written as: 

z' 
x' 

°0 

with x' of weight d. And C' 2 = C' Q + z'. Let c' be in C ' and 
denote by c the extension of c' with (00), then c G Co- Now 
y' ■ c' = since y ■ c = and y' • x' = since y ■ x = 1, 
which proves that for c e C' Q , c ■ y' = 0. Moreover since 
y • z = 0, we deduce that y' ■ z' = 1. The latter results show 
that C is a [An, 2n] self-dual code with minimum distance d 
(since x' has weight d), such that C — C' Q U (z' + C' ) and 
with shadow S' = (y' + C' a ) U (y' + z' + C' ). Finally we 
remark that by construction, for any vector of S' it is possible 
to add either (11), (01), (10), or (00) such that the extended 
vector is in S. Since all the weights of 5" are congruent to 
s — 1 (mod 4) and since y' of weight s — 1 is in S' we 
deduce that the minimum weight of S' is s — 1 which proves 
the theorem. □ 



B. Classification of s-extremal [38,19,8] codes 

Let C be an extremal self-dual [38, 19, 8] code. If C satisfies 
W\ in equation ©, then S\ in (O is also satisfied. So we have 
d = 8 and s = 7; hence 2d+ s = 23 = n/2 + 4. This implies 
that C is an s-extremal code with parameters (7,8). Clearly 
if C satisfies W2, then C cannot be an s-extremal code since 
2d + s — 19 < 23. The s-extremal code can be obtained 
directly from the classification of all [38, 19, 8] by a simple 
computation on the weight enumerator. We obtain: 

Theorem 9: There are exactly 1730 s-extremal [38,19,8] 
codes. 



C. Classification of s-extremal [38, 19, 6] codes 

The case of d — 6 was mainly considered in |2|, where s- 
extremal codes are known to exist for the lengths 22 < n < 44. 
Two lengths 36 and 38 remained open in 0. Later, s-extremal 
codes of length 36 and d = 6 were classified in [1|. The only 
open case is the classification of s-extremal codes of length 
38 and d — 6. There are at least two such codes as shown 
in 0. We show that there are exactly two s-extremal codes 
of length 38 and d = 6. 

For a self-dual [38, 19, 6] code to be s-extremal, the min- 
imum weight of its shadow must be s = 11. A simple 
approach to find all s-extremal [38, 19, 6] codes is to apply 
the recursive construction, starting from the set of all inequiv- 
alent [36, 18, 4] self-dual codes. Unfortunately, since there are 
436, 633 [36, 18, 4] self-dual codes, such a computation would 
require more than 80 days, and although it is doable theoreti- 
cally, in practice it remains largely too costly. Fortunately, by 
using the fact that such an s-extremal code has a shadow with 
high minimum weight it is possible to dramatically decrease 
this computation. 

We have shown in Theorem [8] that it is possible to relate 
the weight of a shadow of a code to the that of the shadow 



of the subtracted code under certain conditions. We use this 
result to prove the following classification theorem: 

Theorem 10: There are exactly two s-extremal [38, 19, 6] 
codes. 

Proof. Let C be an s-extremal [38, 19, 6] code, then C has 
shadow weight s = 11. Applying Theorem [8] we deduce 
that there exist two coordinates on which the subtraction of 
(11) of C produces a [36, 18,4] self-dual code with shadow 
weight 10. Hence if one applies the recursive algorithm 
starting from the set 6*36,10 of all inequivalent [36, 18, 4] 
self-dual codes with shadow weight 10, we construct the set 
of all [38, 19, 6] self-dual codes (up to permutation) which 
by a subtraction of (11) on certain two columns give the 
set 5*36, io- Hence applying the recursive algorithm to S^io 
gives a set of self-dual codes which contains all s-extremal 
[38, 19, 6] codes. In practice, from the classification of ifPTIl . 
there are exactly 24 [36,18,4] self-dual codes with shadow 
weight 10. The application of the recursive algorithm is then 
fast with these codes and we have that there are exactly two 
s-extremal [38, 19, 6] codes. □ 

The two s-extremal codes 638,1, £38,2 have covering radius 
11 and their generator matrices G(C38,i), G(C38,2) are as 
follows: 

10000000000000000000000011111100001110 
01000000000000000000000011111100000001 
00100000000000000001010101010101111001 
00010000000000000001010101010101000101 
00001000000000000000000101010101110011 
00000100000000000000000101010110111111 
00000010000001000000000001111000000000 
00000001000001000000000001110111000000 
00000000100001010000000110100011000000 

<3(C 3S !) = 00000000010001010000000110010000000000 
00000000001001010000010000100101010110 
00000000000101010000010011011010010110 
00000000000011000000000011000011000000 
00000000000000110000000011001100000000 
00000000000000001001000101010110011001 
00000000000000000101000110101001101001 
00000000000000000011000011111111001100 
00000000000000000000110011111111111100 
00000000000000000000001111001111000000 



G(C 38 , 2 ) = 



10000000000000000100010101010110010010 
01000000000000000100010101010110011101 
00100000000000000101000000000000100110 
00010000000000000101000000000000011001 
00001000000000000000000101010101000000 
00000100000000000000000101010110111100 
00000010000000010101010001011101110101 
00000001000000010101010001010010110101 
00000000100000000101010110000110110101 
00000000010000000101010110110101110101 
00000000001000000000000110101001110101 
00000000000100000000000101010110000110 
00000000000010010101010011100110110101 
00000000000001010101010000100101110101 
00000000000000110000000011001100000000 
00000000000000001100000011111111110000 
00000000000000000011000011111111111111 
00000000000000000000110011111111001111 
00000000000000000000001111001111000000 



Notice that these codes were already known from [2|, but 
it was not known whether there exist other codes. 

D. Up-to-date tables for s-extremal codes 

In the following we give up-to-date tables for s-extremal 
codes of minimum distance 6 and 8: 

• d = 6 

For this minimum distance, we know that there are exactly two 
s-extremal codes of length 38 and d = 6 from Theorem [10] 
This was the only unknown case (see Q], 0)- Now we 
complete the classification of s-extremal codes of d — 6 in 
Table HV] 

• d = 8 

In this case, s-extremal codes exist for 32 < n < 44. More 
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TABLE IV 

Number of s-extremal codes with d = 6 



n 


num 


n 


num 


22 


1 


34 


17 


24 


1 


36 


5 


26 


1 


38 


2 


28 


2 


40 


1 


30 


9 


42 


1 


32 


19 


44 


1 



precisely, s-extremal codes of length 32 were known from the 
classification of extremal self-dual codes of length 32, and 
s-extremal codes of length 36 were done in QJ. We have 
completed the classification s-extremal codes of length 38 and 
d = 8 from Theorem [9] We list currently known codes for 
d = 8 in Table [Vj 

TABLE V 

Number of s-extremal codes with d = 8 



n 


num 


ref 


32 


3 


121 


36 


25 


ffl 


38 


1730 


this paper 


40 


> 4 


|9|,|7 


42 


> 17 


[91,16-1 


44 


> 1 


ED 



VI. Covering radii of self-dual codes of length 38 

The covering radius p(C) of a code C is the smallest integer 
R such that spheres of radius R around codewords cover FJ. 

The following theorems give the lower and upper bound of 
p(C) for a self-dual code over GF(2). 

Theorem 11: ((8), Theorem 1) Let C be a self-dual code 
of length n over GF{2). Then n Ci > 2 n ' 2 . More pre- 

cisely, 



2i< P (c) n C 2i > 2 (n/2) 1 and X^+i<p(c) n<?2i+i > 

2 (n/2)-l > where nC . 

means n choose i 

Theorem 12: fll8l.ll20l. Delsarte's bound) Let C be a self- 
dual code of length n over GF(2) and s be the number of 
distinct nonzero weights in C. Then p(C) < s. 

By Theorem [TT| any self-dual [38, 19] code has covering 
radius at least 6. On the other hand, the weight enumerators (0]i 
and © of any self-dual [38, 19, 8] code has 13 nonzero 
weights. Thus by Theorem [12] the covering radius of any self- 
dual [38, 19, 8] code is at most 13. Combining both, we have 

6 < p(C) < 13 for any self-dual [38, 19, 8] code C. 

Using our classification of all self-dual [38, 19, 8] codes, we 
have the following. 

Theorem 13: All 2744 self-dual [38, 19,8] codes have cov- 
ering radius 7. 

Remark 14: If we choose a coset representative of weight 

7 and using it as a vector x in Theorem [2] then the built 
code will be an extremal self-dual [40,20,8] code. Hence for 
n = 38, any extremal self-dual [38, 19, 8] code can produce 
an extremal self-dual [40,20,8] code using the building-up 
construction. This is not always true for some lengths (e.g. 
n = 24, n = 32). 



Proposition 15: Let C\ be a self-dual code of length n and 
covering radius p(Ci). Then any self-dual code C2 of length 
n+2 obtained by the building-up construction (in particular, by 
the recursive algorithm) has covering radius p(C2) < p(C\) + 
2. 

Proof. Let p{C x ) = r. We recall [20, Theorem 1.25.5] that the 
covering radius p{C\) of a linear code C\ with parity check 
matrix H\ is the smallest number s such that every nonzero 
syndrome is a combination of s or fewer columns of Hi, and 
some syndrome requires s columns. The generator matrix G2 
of C2 by the building-up construction is of the form d2j. This 
G2 is also a parity check matrix of C2 as C2 is self-dual. Any 
syndrome u = [u\,U2,--- ,«( n /2)+i] T with respect to G2 
can be written as ui[l,0, ••• , 0] + [0, u 2 ,--- , u^ n / 2 )+i} T ■ 
Now [0,M2,-- - ,U( n /2)+i] T is a linear combination of r 
or fewer columns of G2 as G\ has covering radius r, and 
[1,0, ••• ,0] T is the difference of the first columns of G2 
in the form (0. Hence u is a linear combination of at most 
r + 2 columns of G 2 . Thus p(C 2 ) < p{C\) + 2. □ 

Using Proposition Q3] we have a better upper bound for the 
covering radius of a self-dual [38, 19, 6] code than Delsarte's 
bound as follows. 

Corollary 16: The covering radius p(C) of any self-dual 
[38,19,6] code is 6 < p{C) < 12. 

Proof. The lower bound is true for any even [38, 19] code by 
Theorem QT] Delsarte's bound would imply p(C) < 15. For 
a better upper bound, we recall that any self-dual [38, 19, 6] 
code can be constructed from a self-dual [36, 18,4] code by 
the recursive algorithm. Since the covering radius of any 
self-dual [36, 18, d = 4, 6] code is at most 10 E), it follows 
that the covering radius of any self-dual [38, 19, 6] code is at 
most 12 by Proposition [TBI □ 
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APPENDIX 

Let i = 6,9,12,14,18,21,24,36,144,168,216,342,504. 
Then C7(C| 8 ) represents a generator matrix of a new self- 
dual [38, 19, 8] code C\ & with the automorphism group order 



|Aut(C : 



10000000000000000001010110011101101001 
01000000000000000001010110011101010110 
00100000000000000000010101101000011000 
00010000000000000000010101100111100111 
00001000000001000001100111011001010100 
00000100000001000001100111100110101000 
00000010000001000001110000100101100001 
00000001000001000001110011011010101101 
00000000100001000000100101011001101110 
00000000010001000000100110011010100001 
00000000001001000001000110010000010010 
00000000000101000001000100111100001100 
00000000000011000000000010101111100001 
00000000000000100000110000010000110101 
00000000000000010000000011010011110110 
00000000000000001001100010001100101000 
00000000000000000101010001111111011011 
00000000000000000011110010011111100010 
00000000000000000000001101011111101101 



10000000000000000000011101 
01000000000000000000011101 
00100000000000000001010011 
00010000000000000001010011 
00001000000000000001100001 
00000100000000000001100001 
00000010000000000000111101 
00000001000000000000111110 
00000000100000000001101000 
00000000010000000001101011 
00000000001000000000001011 
00000000000100000000000100 
00000000000010000001001111 
00000000000001000001000000 
00000000000000100001111011 
00000000000000010001001000 
00000000000000001000100100 
00000000000000000100010111 
00000000000000000011111111 



010111111101 
010111000010 
111101100010 
110010010001 
101110010111 
010001100111 
011000110110 
100111111001 
100100111001 
100111111010 
101101000101 
011110110101 
010010110110 
100010110110 
011010100010 
011001101110 
011001010010 
101010101101 
000000001100 



G(Cit) = 



10000000000000000000110110010111100110 
01000000000000000000110110010111011001 
00100000000000000000101110101111010001 
00010000000000000000101110101100101101 
00001000000001000000011101100111001100 
00000100000001000000011101101000000000 
00000010000001000000010000010011110011 
00000001000001000000010000101111001111 
00000000100000000000000111111001011110 
00000000010000000000000100110110010001 
00000000001001000000010001110101011110 
00000000000101000000011101110110010010 
00000000000011000000001100110000111100 
00000000000000100000111000101011010111 
00000000000000010000001011010111100111 
00000000000000001000111101010001111001 
00000000000000000100001110101101111010 
00000000000000000010100101111001101101 
00000000000000000001011001001010101101 



10000000000000000011000110001010110100 
01000000000000000011000110001010001011 
00100000000000000011010110100010000101 
00010000000000000011010110100001001010 
00001000000000000100000111000001101101 
00000100000000000100000111001101100010 
00000010000000000010010000100101101011 
00000001000000000010010011011010010111 
00000000100000000111000011111010101011 
00000000010000000111000100101100111010 
00000000001000000001000101001010011110 
00000000000100000001010001110011111010 
00000000000010000000000110011001011000 
00000000000001000011000110100101010100 
00000000000000100101000101010101010111 
00000000000000010110010010010011001111 
00000000000000001111010000011100000110 
00000000000000000000100000101100110110 
00000000000000000000001000101010101101 



lOOOOOOOOOOOOOOOOOOl 11010001 1110101001 
01000000000000000001110100011110010110 
00100000000000000001010110110111011011 
00010000000000000001010110111000011000 
00001000000000000001000101010001101101 
00000100000000000001000101101110100010 
00000010000000000000011100101101111010 
00000001000000000000011111010010000110 
00000000100000000001001001010001111010 
00000000010000000001001010010010110110 
00000000001000000000101010011000001010 
00000000000100000000100101101011110101 
00000000000010000001101110100111111001 
00000000000001000001100001010111000101 
00000000000000100001111110010000101011 
00000000000000010001001101010011010100 
00000000000000001000100001010011101011 
00000000000000000100010010100000010111 
00000000000000000011111111000000111100 



10000000000000000000101001101100111000 
01000000000000000000101001101100000111 
00100000000000000000111010000100010100 
00010000000000000000111010000111100111 
00001000000000000000010111111000111111 
00000100000000000000010111110111001111 
00000010000000000000110101010011000000 
00000001000000000000110101100011111100 
00000000100000000001111110100110000110 
00000000010000000001111101010110110101 
00000000001000000000001010111110011011 
00000000000100000000000110000010011011 
00000000000010000001001110000010010100 
00000000000001000001000001001101100111 
00000000000000100001010010011000100010 
00000000000000010001100010010111101101 
00000000000000001000001110011011011110 
00000000000000000100111101100100100001 
00000000000000000011111100110000001111 



lOOOOOOOOOOOOOOOOOOOl 101010001101111 10 
01000000000000000000110101000110000001 
00100000000000000000010101101110110011 
00010000000000000000010101101101111100 
00001000000000000000010111111011001100 
00000100000000000000010111110111000011 
00000010000000000000110101010011110011 
00000001000000000000110101100000001100 
00000000100000000001101100100101101011 
00000000010000000001101111010101010100 
00000000001000000001110001100101101110 
00000000000100000001111101011001010010 
00000000000010000000110101011001101110 
00000000000001000000111010010101011101 
00000000000000100000111011000011000110 
00000000000000010000001011001100001010 
00000000000000001001100111000011001010 
00000000000000000101010100111100110101 
()()()()()()()()()()()()()()()()()()1111110011()()11 111100 



10000000000000000100000111100100111000 
01000000000000000100000111100100000111 
00100000000000000101000100010111001111 
00010000000000000101000100011000110011 
00001000000000000000100101101000100001 
00000100000000000000100101010111010001 
00000010000000000101100000010010111010 
00000001000000000101100011101101111001 
00000000100000000100110101101110001001 
00000000010000000100110110101101001010 
00000000001000000101010110100111111001 
00000000000100000101010111110010100010 
00000000000010000100010010011000000101 
00000000000001000100010011001110101101 
00000000000000100101100010010110100100 
00000000000000010101010001010101100100 
00000000000000001100110011110011110000 
00000000000000000011110001100110100100 
00000000000000000000001110100110101011 



10000000000000000001000010111010001000 
01000000000000000001000010111010110111 
00100000000000000101000100010111001111 
00010000000000000101000100011000001100 
00001000000000000001001100001100010100 
00000100000000000001001100110011100111 
00000010000000000100000111010000100111 
00000001000000000100000100101111011011 
00000000100000000101001100111011001010 
00000000010000000101001111111000001001 
00000000001000000100101111110010110101 
00000000000100000100100000000001110101 
00000000000010000101101011001101001001 
00000000000001000101100100111101110110 
00000000000000100101101100110000001111 
00000000000000010101000001100100100010 
00000000000000001100101101100100101110 
00000000000000000011100001010111011110 
00000000000000000000011110010111101101 



10000000000000000100001010111101011100 
01000000000000000100001010111101100011 
00100000000000000100100101000100101101 
00010000000000000100100101001011011110 
00001000000000000100101001000010111001 
00000100000000000100101001111101001001 
00000010000000000000000011001101011011 
00000001000000000000000000110010011011 
00000000100000000000110111100010110101 
00000000010000000000110100100001111010 
00000000001000000000110101111000100111 
00000000000100000000111010001011010100 
00000000000010000000010000010100111001 
00000000000001000000011111100100111001 
00000000000000100100001101100011101110 
00000000000000010100111110100000100010 
00000000000000001100110011110011110000 
00000000000000000010011110010011101110 
00000000000000000001100001010011101101 



10000000000000000000001101101111100110 
01000000000000000000001101101111011001 
00100000000000000000111010000100101101 
00010000000000000000111010000111100010 
00001000000000000010011111111110011011 
00000100000000000010011111110010010100 
00000010000000000010111000100011001111 
00000001000000000010111011011100111100 
00000000100000000000001000111001010001 
00000000010000000000000111000101101101 
00000000001000000000111110110000110110 
00000000000100000000000010001100111010 
00000000000010000000101110011111110011 
00000000000001000010011101110000111010 
00000000000000100010001110010110100010 
00000000000000010000000010000110010111 
00000000000000001000010011110110100010 
00000000000000000110010000010101100100 
OOOOOOOOOOOOOOOOOOOl 11001 1010011111001 
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10000000000001000000001010101111111011 
01000000000001000000001010101111000100 
00100000000000000100001101101001111010 
00010000000000000100001101100110111001 
00001000000001000100101001000000101011 
00000100000001000100101001111111100111 
00000010000001000000101011100010011110 
00000001000001000000101000011101101110 
00000000100001000000101000000110110011 
00000000010001000000101011000101111100 
00000000001001000000000010110001111010 
00000000000101000000001101000010000110 
00000000000011000000001111110000110000 
00000000000000100100100101001110000101 
00000000000000010100001001101011101101 
00000000000000001100101100010101100111 
00000000000000000010101001011000010001 
00000000000000000001001001111110110101 
00000000000000000000011111100110010111 



[20] W.C. Huffman and V. Pless, Fundamentals of Error-Correcting Codes, 

Cambridge University Press, 2003. 
[21] J.-L. Kim, "New extremal self-dual codes of lengths 36,38 and 58, IEEE 

Trans. Inform. Theory, vol. 47, no. 4, pp. 1575-1580, 2001. 
[22] F.J. MacWilliams, N.J. A. Sloane and J.G. Thompson, "Good self-dual 

codes exist," Discrete Math. vol. 3, pp. 153-162, 1972. 
[23] G. Nebe, E. Rains and N.J. A. Sloane, Slef-Dual Codes and Invariant 

Theory, Springer, Berlin, 2006. 
[24] V. Pless, "A classification of self-orthogonal codes over GF(2)," Dis- 
crete Math. vol. 3, pp. 209-246, 1972. 
[25] E. M. Rains and N. J. A. Sloane, " Self-dual codes," in Handbook of 

Coding Theory, ed. V. S. Pless and W. C. Huffman. Amsterdam: Elsevier, 

pp. 177-294, 1998. 



10000000000000010010011000100100010110 
01000000000000010010011000100100101001 
00100000000000010010010110111100011000 
00010000000000010010010110110011010100 
00001000000000000010011001010011110000 
00000100000000000010011001101100001100 
00000010000000010010000100000110000011 
00000001000000010010000111111010110000 
00000000100000000000010111100110101101 
00000000010000000000011000100101101110 
00000000001000010000010100011110100100 
00000000000100010000011111010101011101 
00000000000010010000010010001010111100 
00000000000001010010000000010101100010 
00000000000000110010011010100111100111 
00000000000000001000001100010110101110 
00000000000000000110011001111110111001 
00000000000000000001010010100011010010 
00000000000000000000110111001000000101 
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